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I. INTRODUCTION 



Spacetime fluctuations are expected to be a generic feature of a theory which combines 
quantum theory and gravitation. While there is still no complete quantum theory of gravity, 
it is possible to investigate some of the characteristics expected of fluctuating spacetime 
geometries. One can roughly classify these fluctuations as being either active or passive. 
Active fluctuations arise from fluctuations of the dynamical degrees of freedom of gravity 
itself, that is, from the quantization of gravity. Passive fluctuations arise from fluctuations 
in the stress tensor due to quantum matter fields. In general, one expects both types of 
fluctuations to be present. One approach to the study of fluctuating spacetimes is stochastic 
gravity pj. More generally, there has been considerable activity in recent years in the 
area of quantum gravity phenomenology, which seeks to find observational signatures of 
the quantum nature of spacetime |E 0, IE 0-0-0-0- E3- EE 12> Jll|- There has also been 
considerable attention given to the effects of classical stochastic gravitational fields k! 3 



lR IE 13 an d to scattering of probe particles by gravitons in an S-matrix approach lj. 

Spacetime geometry fluctuations should in principle produce observable effects on test 
particles, such as light rays. Several effects of fluctuating gravitational fields on light prop- 
agation have been discussed by previous authors. For example, Sachs and Wolfe 1*2 



treated the scattering of cosmic microwave photons from cosmological density perturba- 
tions. Zipoy argued that there will be apparent luminosity variations in a source seen 
through a stochastic gravitational field. More recently, this effect has been treated [l2| us- 
ing a Langevin form of the Raychaudhuri equation |2j. In the latter approach, luminosity 
fluctuations are a signature of passive spacetime geometry fluctuations caused by quantum 
stress tensor fluctuations. Zipoy |3] also examined changes in angular position due to a 
stochastic gravitational field. Kaufman [15] considered redshift fluctuations from a bath 
of gravity waves. If one goes beyond the geometric optics approximation to consider elec- 
tromagnetic wav es p ropagating in a stochastic gravitational field, there will be phase shift 
variations [E EE 13 • 

The purpose of the present paper is to analyze in more detail two particular signatures 
of spacetime geometry fluctuations, redshift fluctuations and angular blurring of images. 
We derive simple expressions involving the Riemann tensor correlation function and use 
these to examine fluctuations in redshift and angular position of a source. We will illustrate 
this approach using the cases of gravitons in a squeezed state and in a thermal state, and 
calculate the broadening of spectral lines and an angular blurring of an object viewed through 
a region of spacetime filled with gravitons. This will serve as a simple example of spacetime 
undergoing active quantum fluctuations. 

This paper is organized as follows: In Sect. |H] we provide a detailed derivation of the 
equations for redshift fluctuations and angular blurring given in Ref. [2(j. The derivation 
is independent of the source of fluctuations, being based on the Riemann tensor correlation 
function. In Sect. Illll we assume that gravitons in a general squeezed state act as the source 
of fluctuations and are described by the correlation function, and we then specialize the 
results from Sect. |H]to this model. Several examples are calculated for the variance of both 
redshift and angular position, and are compared with the expected classical time variation of 
these quantities. We then express these results in terms of the energy density of the gravity 
wave and make an order of magnitude estimate for the size of the effect. In Sect. IIVI we 
assume a thermal bath of gravitons as the source of fluctuations and again specialize the 
results of Section [TT1 to this model. The results of the paper are summarized and discussed 
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in Sect. |V] A brief review of the pertinent information regarding squeezed quantum states 
required for this paper can be found in the appendix. Units in which G = h = c = 1 will be 
used, where G is Newton's constant. 



II. REDSHIFT AND ANGULAR BLURRING IN LINEARIZED GRAVITY 



A. Linearized Redshift 

Let be the 4-velocity of a source and v M the 4-velocity of a detector at the events of 
emission and absorption, respectively, and let k 11 be tangent to the null geodesic connecting 
them (see Fig. ^ e.g., path DA). Let A be an affine parameter that runs from A = 
(emission) to A = A (absorption). Since /c M is tangent to a geodesic, 

^ + = o (i) 

and so 

P(A) = F(0) - / d\> r£ (\')k a (\')kP(\'). (2) 
Jo 

This is an integral equation for fc M (A), but in the linearized theory, we let k^{X') ~ fc M (0) on 
the right hand side. So at the point of detection, 

fc"(A ) = fc"(0) - [ X ° d\'T^(X')k a (0)k^0). (3) 



JO 

The frequency of emission is proportional to —k^t^. The proportionality constant depends 
on the choice of affine parameter. If we choose the affine parameter such that, for example, 
k^ = (u, a), 0, 0) in the rest frame of the source, then uj = —t fJj k tJl . Since Eq. (JTJ) is scale 
invariant, we may also select an affine parametrization such that k^ = (1,1,0,0). This 
amounts to choosing the affine parameter to coincide with the source's proper time at the 
evaluation point. We adopt the latter parametrization of A, in which case the ratio of the 
detected to the emitted frequency is 

= -^P(Ao) = -v^(0) + v, f X ° dX' TUX')k a (0)k^0) , (4) 
^0 Jo 

where Uq = <^(0) is the emitted frequency. The first term on the right hand side of this 
equation can be written as 

_^P(0) = -^P(O) - - t M )P(0) = 1 - K - ^)P(O). (5) 

With this, the fractional change in detected frequency between source and detector is 



— = ^ (Ao) U ° = -K-WI+b, / °dAT^(A')r (0)^(0). (6) 



The first term on the right hand side is a Doppler shift which may not be zero even in 
flat space. The second term is a linearized gravitational redshift which depends on the 
intermediate geometry between source and detector. We will focus our attention on this 
second term. 



3 



A 




FIG. 1: A source moves along a worldline with tangent while a detector a proper distance s away 
moves along a worldline with tangent . The source emits a ray at point D which has tangent 
k^(X = 0) at point D and tangent k^(Xo) at A. Parallel propagation of k^ around ABCD results 
in a slightly rotated vector k^ + Ak^. The closed path ABCD encloses the spacetime region of 
interest. 

B. Rate of Change of Redshift 

The next step is to examine the rate of change of the redshift. We examine the behavior 
of successive photons emitted from the source. Our interest lies in the change in frequency 
due to the effects of gravity rather than as a result of a variation of the output of the source. 
This condition is enforced by requiring that the values of k 11 at the start of any two null 
geodesies are related by parallel transport along the worldline of the source. Thus we assume 
that Uq is constant between successive emissions and want to find cu(X ) at the detector due 
to changes in spacetime geometry. Initially at the detector, 



is the frequency at the detector at proper time t\. While after a time At, at proper time 
r 2 , the frequency at the detector is 



w(ti, Aq) 



(7) 



ufa, Aq) 



v^r 2 )k^{T 2 , X )u . 



(8) 



To find v^(r 2 ), we first note that t> M is tangent to a geodesic: 



dv» 



(9) 



dr 



+ 



so 



(r 2 ) - W "(7i) 



dv^ 



At ~ -r^^/Ar. 



(10) 



v 
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This depends on F^g at the location of the detector. For simplicity, we assume F^g = at 
the location of the detector. This can be achieved by assuming the detector is located in a 
flat region; more will be said about this in the next subsection. Thus v^(t2) = v^{t\) = -iA 
To find the change in we parallel transport the vector around the closed path ABCD 
(see Fig. |TJ). Recall that if we parallel transport a vector around a closed path, the change 
in the vector can be expressed as an integral of the Riemann tensor over the area enclosed 
by the path. For parallel transport around an infinitesimal parallelogram, we have for an 
arbitrary vector V 1 

V»^V> t + 6V» , (11) 

where 

5V» = -R" avp V a t v k p AXAT. (12) 
Integrating over A and r transports V 1 around a finite parallelogram so that 

AV = - j dr J dXR\ u(3 V a t v k p . (13) 

Since we assume the detector is located in a flat region, we can further assume that k^ may 
be trivially transported from point A to point B, and we have 

F(r 1 ,Ao)+A^ = F(r 2 ,A ). (14) 

Thus ^ 

A^ = ^(r 2 ,A )-^(r 1 ,A ) = - I * dr I ° ' dXR^k^k^ , (15) 

J-n Jo 

and from Eqs. (J7J) and (JSJ) 

A^(Ao) oj(r 2 ,Xo)-oj(r u X ) , 1U f T2 , f x ° 



U UJ 



-VfAk" =v lt ( 2 dr [ ° dX R\ up k a t v k p . (16) 

Jt, jo 



This equation can be thought of as relating the difference between taking path DAB and 
path DCB. Taking limits in the previous equation yields the rate of change of the fractional 
redshift in the linearized theory: 

^Hf dAfl -^ s (17) 

An equivalent formula has been obtained by Braginsky and Menskii [16 : ] using the geodesic 
deviation equation. 



C. Fluctuating Redshift 

Now suppose the Riemann tensor is subject to fluctuations - active, passive, or both. We 
first have to specify and and how they behave under the fluctuations. The simplest 
assumption is that f 1 and do not fluctuate, and in the underlying flat geometry t> M = t M . 
Physically, this might be achieved via several methods. We have already mentioned that we 
can assume the detector to be located in a flat region. However, one could also consider the 
source and the detector rigidly attached to one another by non-gravitational forces. The 
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same effect might be achieved if the source and detector are separately attached to platforms 
(e.g. planets) which are large enough that they travel on an average geodesic in the mean 
spacetime. This can happen if the spatial average of the fluctuations over the platform is 
small. For our purposes it is sufficient to assume the perturbation vanishes at both the 
source and detector, e.g., a gravity wave passes between source and detector, but far from 
either source or detector; so we can assume the source and detector to be located in flat 
regions. The fluctuations of the Riemann tensor are described by the correlation function 

CaPuviS/xrfax) = (Ra/3lJ,u(x)R 7 Spa(x )) ~ {R a /3fiu(x)) (RjSpa(x')} , (18) 

where the indices <xfi\w refer to point x while the indices jdpa refer to point x'. With 
v 1 * = t M , we get from Eq. (J§J) 

e - — = = f Xo dX K(x)k * (o^o) v (19) 

^0 ^0 Jo 

We let r^g fluctuate with fixed k a = k a (0) and t M , so that 

(0= [ X °d\(T^(\))k a k\. (20) 
Jo 



The integral in this equation is evaluated along a single line, e.g., along the paths DA or CB. 
The equation may be interpreted in the following way: Suppose the spacetime is subject to 
quantum fluctuations. Then, given an ensemble of systems, measurement of the fractional 
redshift along the same line will yield different results. This equation gives the expectation 
value of those measurements. Therefore, comparing the result of an integration along path 
DA with that along path CB will give the expectation of the difference of the fractional 
redshifts. Another, more convenient way to obtain the same information is via Eq. 1)1 7]). 
Notice that 

d£ d fuj(X )-uj \ d f uj(X )\ ^ 



dr dr \ u J dr \ u 

and thus from (JTTJ) , 

di 
dr 

Integrating this expression yields 



I" d\R a ^{T,\)t a kH»k v . (22) 
Jo 



[ 2 dr [ ° d\ R al3 pv{T, X)f*k p lPkr. (23) 
Jti Jo 



We now let Rap^vij, A) fluctuate, and find the expectation value (A£) is 

(Af) = f 2 dr 1° d\ (R a ^{r, X))t a kH»V . (24) 

Jti Jo 

The expectation of the square, ((A£) 2 ), is related to the line broadening an observer 
will see. Recall that £ is the fractional change in frequency of an observed spectral line, 
the fractional redshift. Upon quantizing the perturbation between source and detector we 
strictly must discuss the fractional redshift along a given path as an ensemble average. An 
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observer, over some proper time interval At, collects information on a distribution of frac- 
tional redshifts; ((AO 2 ) is the squared width of this distribution. The physical realization 
of this measurement is a broadening of the observed spectral line. 

However, there is a contribution to spectral line broadening due to regular time dependent 
variations of the spacetime, for example as a result of passing classical gravity waves. We 
thus characterize the fluctuation of redshift about the classical time dependent variation 
that arises from a nonzero expectation value of the Riemann tensor. Using Eq. (|18|). we can 
express the variance of the fractional redshift, 5£ 2 , as 

5£ 2 = ((AO 2 ) - (AO 2 = J da J dd ' C a ^ vl&pa {x,x')t a kH p k v Vkh p k° ". (25) 

The integration over the spacetime region J da corresponds to dr J Q A ° dX and similarly 
for da'. This is an integration over the spacetime region enclosed by the path ABCD in 
Fig. The point x corresponds to the point (r, A) and similarly for x'. 



D. Fluctuating Angular Position 

We can also relate the degree of angular blurring of the source observed by the detector 
to the Riemann tensor correlation function. Let s M be a unit spacelike vector in a direction 
orthogonal to the direction of propagation of the null rays; thus = s p k p (X = 0) = 0. 
Then at the observation point 

s^{X ) = tan 6 « 9, (26) 

where G is an angle in the plane defined by the pair of spacelike vectors and n p = k^ — t p 
and is assumed to be small, |G| <C 1. The angle G is the angular deviation in the direction 
of of the image of the source from its classical flat space position. A treatment similar to 
that shown for fractional redshift allows us to express the change in angle, AG, in terms of 
an integral of the Riemann tensor as 

AG = S(1 AF = J daR a ^ v s a k^k v . (27) 

A fluctuating spacetime results in an ensemble distribution of image positions about the 
classical flat space position. Analogously to the line broadening effects, the fluctuating 
angular position manifests itself as a blurring of the source's image. ((AG) 2 ) is therefore a 
measure of the angular size of the image. The variance of AG, <5Q 2 , due to fluctuations in 
the Riemann tensor is 

5Q 2 = ((AG) 2 ) - (AG) 2 = J da J dd C a ^ ul&p(7 {x,x')s a kh p k h 's<k s t p k° . (28) 



E. Quantization of the Riemann Tensor 

Our key results for the redshift fluctuations, Eq. (|25|). and for the angular blurring, 
Eq. (J2EJ), apply both to active and passive fluctuations of spacetime geometry. We will 
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present some explicit examples for the case of active fluctuations. We examine fluctuations 
produced by gravitons occupying a squeezed state and then a thermal bath of gravitons in 
the following two sections, respectively. First we describe some of the formalism relevant to 
both cases. A linearized quantum field theory for gravity is used, with the field operator 
expanded as 

V = J2 ( A ^P eUrXT + h - c ) > ( 29 ) 
e,p 

where £^ is the wavevector and p labels the polarization of a mode with polarization tensor 
A^y = A^ijp). In this theory, the Riemann tensor operator is given by 

R a p^u{x) = d v d[ a hp]fj, — dfj,d[ a hp] v . (30) 



Here the convention for antisymmetrization as found in Ref. |21[ is used, i.e. d v d[ a hp\ p = 
\{d u d a hp p — d v dph ap ). Henceforth and R a pp,u are understood to be operators, and the 
hat may be suppressed without confusion. The expectation value of R a pp, u is 

(Raf3fiv{x)) = (d u d[ a h^) - (dpdlahflv) . (31) 

It is convenient to define 

K af}p , vlSp(T (x,x') = (d u d a h f3fl (x)d' (T d^hs p (x')} - {d v d a h^{x))(d' a & 1 hs p {x')) 

= dudaff^ {(hfr(x)h Sp {a/)) - (hfr(x))(h Sp (a/))) , (32) 

where d' denotes differentiation with respect to x'. The Riemann tensor correlation function 
may be expressed as 

Capfiv-ySpaix, x) = 4^ Q/ g] [^] [ 7 S][ p(T ] (x, x) . (33) 



III. GRAVITONS IN A SQUEEZED STATE 



A. Single Mode Squeezed State 

Gravitons in a squeezed state are the natural result of quantum graviton creation in a 
background gravitational field, such as in an expanding universe [22] . Here we suppose that 
the region between a source and a detector is filled with gravitons in a squeezed state, which 
produce spacetime geometry fluctuations. A brief summary of the properties of squeezed 
states needed for the following calculations may be found in the appendix. To calculate 
Caf3pu-y8pa(%, %'), we must calculate {hp^,{x)hs p {x')) and (h/3 p (x)). We are interested in the 
change in quantum fluctuations between the vacuum state and a squeezed state, so we may 
take the correlation function to be normal-ordered. We choose to evaluate the expectation 
value with respect to a gravity wave in a general single mode squeezed state \a,() (see the 
appendix), where a and ( are the displacement and squeeze parameters, respectively, and 
the only mode excited has a specific wave vector frequency u g = £° and polarization 
p. We may now proceed to calculate ((,a\h pu \a,(), (C^ a \hfj,u(x)h a /3(x')\a, () and then find 
: K a ^ ulSpa (x,x') :. The result is: 

• ^aPfjiu -ySpa (•£; X ) . £ a £i/£^£ (T Ap p Ag p f(^X J x ) , (34) 
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v = u*+2s 




FIG. 2: The spacetime region between the source and detector in null coordinates. For each fixed 
u = u* between u = and u = to, one integrates along the line u = u* from v = u* to v = u* + 2s. 



where 



f{x, x') = [cosh(2r) - 1] cos[£ T (x T - x' T )} - sinh(2r) co&[l T {x r + x' T ) + 0]. (35) 



The parameters r and 6 are defined such that £ = re* 6 (see the appendix). Using Eq. (|33j) . 
the normal ordered Riemann tensor correlation function can be expressed as: 



Cap^spa&x') :=4 (£[ a y4 /3 ][ M 4])(£[ 7 A (5 ][ p 4])/(x,a; / ) 



(36) 



Note that the correlation function depends only upon the squeeze parameter (, and not 
upon the coherent state parameter a. 

We now use this normal ordered version of the Riemann tensor correlation function in 
Eqs. (J2SJ) and When performing the da and da' spacetime integrations on f(x,x'), it 

is convenient to use null coordinates, u = t — x, v = t + x (see Fig. |2J. The result is 



F(u g ,s,t ) 



//•to rt nu+2s ru'+2s 

da'f(x,x')= / du du' / dv dv' f(x, x') 

JO JO Ju Ju' 

16{1 - cos[s(4 - cog)]} [1 - cos(u g t )] 

U 2 g (l x -UJg) 2 

x [cosh(2r) — sinh(2r) cos(^ + s£ x — uj g (t + s)) — 1] . 



(37) 



Here s is the spatial separation of source and detector, and to is the observation time. It 
should be noted that, since we integrate over a spacetime slice of constant y and z, we may 
set y — y' = z = z' = 0. Finally, from Eq. (J2SJ), the redshift fluctuations become: 

(38) 



5e = 4(£ [a ^ ][ ^ ] )(£ [7 A 5][ / CT] )t Q A;^A;^A;VA; CT F(^, S ,t ). 
While from Eq. (|28|) . the angular fluctuations are given by 

5Q 2 =4(f [a A ffl ^])(f [7 i % 4])s Q ^rMVF(a; 9)S ,to). 



(39) 



Interestingly, since F(u g , s, t ) is independent of the displacement parameter, a, so are S^ 2 
and 5Q 2 . Therefore, the fluctuations depend only on the squeezing parameter, (, and we 
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can immediately say that a coherent state (classical wave) for which ( = induces no 
fluctuations. For the following calculations, we assume the passing gravity waves are in 
the Transverse Tracefree (TT) Gauge. However, note that since our results derive from 
the Riemann tensor, the equations are gauge invariant, and we use the TT gauge only for 
calculational convenience. 

To begin, note that if the gravity waves and photons travel in parallel, then one finds 
(£ [a A l3][f ,£ u] )t a k f3 t' M k' / = 0. This indicates that there is no induced redshift fluctuation nor 
angular blurring due to a squeezed state gravity wave traveling with the photons. This 
is as one might expect, since a gravity wave in the TT gauge has only non-zero transverse 
components. We may subsequently restrict our attention to transversely propagating gravity 
waves and will assume the gravity waves to propagate with £^ = u g (l, 0, 0, 1) while the 
photons continue to have wave vector k^ = (1, 1, 0, 0). We define 

$(ujg, s, to) = [1 — cos(o; 5 s)] [1 — cos(a; ff to)] [cosh(2r) — sinh(2r) cos(6 ) — uj g {t + s j) — 1] , (40) 

so that for a gravitational wave propagating in the z direction, for which £ x = 0, 

16 

F(uj g ,s,t ) = — $(uj g ,s,to) . (41) 

LOg 

1. Redshift Fluctuations 

With transversely propagating gravity waves, 

{£ [a A^£ u] )t a k^k v = -\co 2 g A + , (42) 

and therefore 

5e = 4Al$(LUg,s,t ). (43) 

There exists a non-zero effect which in this case depends only on the (+) polarization. For 
a gravity wave propagating in the ^-direction, A + = A xx = —A yy . 



2. Angular Blurring 

a. Case 1 Let the photons propagate in the x direction while the gravity waves prop- 
agate in the z direction and probe the y component of angular blurring, s M = (0, 0, 1, 0). We 
find 

(£ [a A m / u] )s a k^k» = - A u 2 g A x (44) 

and 

SQ 2 y = 4 A 2 x $(ujg, s,t ). (45) 

Note that here the y component of blurring depends only on the (x) polarization of the 
gravity waves. Here A x = A xy = A yx . 

b. Case 2 To probe the z component of blurring, set s M = (0, 0, 0, 1). The result is 

6e 2 z = 4Al$(u g ,s,to), (46) 

which is the same as Eq. f!43|) for redshift fluctuations, and only depends on the (+) polar- 
ization of the gravity waves. 
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B. Classical Time Dependent Variation 

As noted above, the fluctuations in redshift and angle depend upon the degree of squeez- 
ing, measured by the parameter (. In this subsection, we will examine the expectation 
value of the change in fractional redshift, (A£) and in angle, (AO). These quantities will 
depend only upon the displacement parameter a, and hence would be the same in the co- 
herent state |a,0) as they are in the squeezed state \a, £). For this reason, we regard these 
quantities as giving the classical time-dependence. Alternative ap pro aches to the classical 
time-dependence of these quantities may be found in Refs. [3, [Til Il8| 

Here it is possible to calculate (A£) as a single integration of the Riemann tensor over 
the spacetime region da via Eq. (|17|). However, for comparison purposes, we calculate (A£) 2 
directly. In this case, 

(A£) 2 = J da J da' {R a p^{x)){R 1 8 P a{x'))t a kh^k u f i kh p k CJ , (47) 



where now 



with 



(R a ^v(x))(R l5pa (x')) = 4^(£ [Q A /3][M 4 ] )(£ [7 A (5][ / CT] )/ , (x,x / ) (48) 



f(x,x') = a 2 e ie ^ xT+x ' T) + {a*fe- l ^ xT+x ' T) + 2\a\ 2 cos[£ T (x T - x' T )]. (49) 
Integration of f'(x, x') over the spacetime region of Fig. ^ yields 

F'(u) g , s,t ) = J da J da' f'(x,x') = 

16{1 - cos[s(4 - ^)]}[1 - cos(u g t )} f ae , 2[lxS ^ g[s+to)] + a ^ [exS ^ g{s+to)] \ 2 (5Q) 



^ 2 (4-^) 2 

Thus the classical time variation of the redshift is characterized by: 

(AO 2 = 4FV s ,s,to)(^[^/ 3 ][M4])(^ 7 ^][p4])t a A; /3 ^WA;Vr. (51) 

For the case of redshift variations with photons and a gravity wave propagating perpendic- 
ularly (Sect. HI1 A 10 . 

(A0 2 = 4A 2 + $'(uj g ,s,to). (52) 
Here $'(u) g , s, to) is defined in a similar way as $(u) g , s, to) when £ x = 0, 

&{u g , s, t ) = [1 - cos(ov)] [1 - cos(c^o)] («e-^ (s+ ' o) + «*e^ (s+ * o)) ) * . (53) 

The classical time variation of angular position yields similar results. Particularly, for per- 
pendicularly propagating photons and gravity waves and probing the y component of blur- 
ring (case 1 of Sect. IIII A 2J) . one finds 

(AQ y ) 2 = 4A 2 x $'(u Jg ,s,t ). (54) 

While probing the z component of blurring (case 2 of Sect. IIII A2j) . the result is 

(A0 2 ) 2 = 4^'K,s,t o ). (55) 
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Note that the function F'(u g , s,t ), and hence $'(u g , s,to), depends on the displacement 
parameter a, but not on the squeeze parameter, r. Therefore the same can be said for 
results Eqs. f|52jh (J53|) . and (|53j). In particular $'(w g , s,to) = for a = 0. This means a 
coherent state (classical wave), for which r = 0,a^0, exhibits regular time variations but 
no fluctuations. Indeed, recall from Eq. (|37|) that 5£ 2 = 5<d 2 = for r = 0. We will exploit 
this fact shortly by considering a state for which a = 0. 

In a time-averaged measurement, the classical time variation will produce line broadening 
and angular blurring, such as do the quantum spacetime fluctuations from a squeezed vacuum 
state. The two effects can, however, be distinguished in principle. If one were to make 
repeated measurements at the same point in the cycle of a gravity wave, the effects of 
classical time dependence would disappear, but those of quantum fluctuations would remain. 



C. Stress Tensor for Squeezed State Gravity Waves 

We would like to obtain an order of magnitude estimate for ((A£) 2 ) and ((A©) 2 ). How- 
ever, the gravitational wave amplitude A^ v is not directly measurable. We therefore would 
like to express the previous results in terms of the energy density of the waves. While the 
stress-energy tensor is not well defined for a gravitational wave, one can define an effective 
stress tensor in the linearized theory (see, e.g., Ref. 0]). Classically, this effective stress 
tensor is defined as 

= ^ (h^,XL) • (56) 

Here the brackets denote a spatial average over several wavelengths and the TT superscript 
indicates we are working in the Transverse Tracefree gauge. We use this expression with the 
expectation value of a gravity wave in a squeezed state, and perform a spatial average on 
the quantity 

(C,a| : hlp^hl^ : \a,() (57) 
This operation results in an effective stress tensor of the form 

T ^ = W- E ^M^M 2 + cosh(2r)] (58) 

OZ7T z — ' 

a,j3 

while 

: T^y : = — | J^V^M 2 + cosh(2r) - 1]. (59) 
OZ7T z — ' 

The vacuum energy term, Tqq c = T Q0 — : T 0Q : is then 



rpvac 

00 ' 327, 



Examining the contribution from the (+) polarization, where A + = A xx = —A yy and A x 



A xy = A yx = 0, one finds 



T oo C = ^ A K ( 61 ) 
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However, since the quantum vacuum energy is also \oj 9 , we can use the vacuum mode to fix 
the normalization. 

VT ^ = ^ A >l V = r 9 ( 62 ) 
for each £, where V is the quantization volume, which leads to 



A+ = \l^r,- ( 63 ) 



CJgV 



For the (x) polarization, one also finds 



A, = ^. (64) 



The previous results, Eqs. ()43|) . (|45|h and (|4*U|l become 

32tt 

hi 



5? = 6B 2 y = 6e 2 z = — 3(u g , s, t ) . (65) 
to „ v 



So far, we have considered a single mode. However, we may also have a situation in which 
many modes are excited. In this case, we would insert a sum on modes on the right hand 
side of the above equation. If the density of states is large, we let 

Z-jkJ**- (66) 

Suppose the distribution of gravity waves in £-space is narrowly peaked about some u g with 
characteristic width Au) g . If Au g is small, then the integrand is essentially constant over 
the region of integration. The result of the integration is just the integrand multiplied by 
a volume, (A£ x )(A£ y )(A£ z ), in £-space, with A£i a bandwidth in the i-direction of £-space. 
Specifically, 

6e = ^jL (A£ x )(A£ y )(AQ F(u g , Mo)- (67) 



D. Estimate of 5£ 2 and <58 2 

Recall that from Eqs. and (|55)l. the classical time variation depends on the 

displacement parameter, a, but not on the squeeze parameter (. For the purpose of obtaining 
order of magnitude estimates of 5£ 2 and SO 2 , suppose a = and r ^> 1; further assume the 
+ polarization, so that A x = 0. Since for a = 0, $'(ujg, s,t ) = 0, it is clear that (A£) 2 = 
and (AO) 2 = 0, from Eqs. (|52~j) and ([55)1 . respectively. Therefore, for the case a = one 
finds 5£ 2 = ((AO 2 ) and 59 2 = ((A6) 2 ). 

Using Eqs. (joTJj) and and integrating over a sharply peaked distribution function in 
£-space, the energy density for large r becomes 

: T 00 : « £^A£ x M y A£ z . (68) 
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From Eqs. ()37|) and (|67|). we have for large r, 



2e 2r 

((AO 2 ) « A^A^Ai,. (69) 



Thus 



((AO 2 ) -64^1^ = 64^1^41, (70) 



where £pi is the Planck length. With ui g = 2Tr/X g , we have 

((AO 2 ) « -4iA 2 : T 00 : or (A^r™ « ^iA ff ^Tl^7. (71) 

7T 

Suppose for example a gravity wave with X g — 1 ly = 10 18 cm and the closure energy density 
: Too : = 10 8 cm" 4 . Then with £ Pl = l(T 33 cm, 

(AOr ms ~ 10- 33 10 18 VnF = 10~ n . (72) 

By comparing Eqs. (jjS)) and (|4l)|). the order of magnitude of (A6) rms will be the same as 
the result for (A£) rms . 

In principle, the effect can be made as large as desired by increasing the squeezing pa- 
rameter r, which increases the energy density of the wave. However, as shown from the 
estimate given, this is a very small effect and is likely to be unobservable in the present day 
universe. This example serves as a useful model for spacetime geometry fluctuations, which 
will have large effects in the early universe and in the vicinity of black holes. 



IV. THERMAL BATH OF GRAVITONS 

In the previous section, the fluctuation effects of a gravitational wave in a single mode 
squeezed state were examined. Another useful example is a thermal bath of gravitons, 
such as might be created by an evaporating black hole. We now consider such a thermal 
bath as the source of spacetime fluctuations. In particular, suppose the spacetime geometry 
fluctuates in such a way that (T^) = (R^^g) = 0, but (R^^pR 1 s pa ) P ^ 0- 111 effect, we 
are ignoring the average spacetime curvature due to the bath of gravitons. In this case the 
Riemann tensor correlation function is 

C a f3fiwy8pa( x ^ x ) = (Ral3fj,u(x) R^Spa(x ))/3 , (73) 

and therefore 

5£ 2 = ((AO 2 ) = [ da [ da' {R aPl , u {x)R l5p(T {x'))pt a kH^k u f t k 5 t p k a (74) 

and 

6& 2 = ((A0) 2 ) = 

where (R a g^ v {x)R 1 spa{x')) g 
tion. 



= J da J da' {R a0llu (x)R 1&pa (x')) (5 s a kH^k v s 1 k 5 t p k CT , (75) 
is the thermal normal-ordered Riemann tensor two point func- 



14 



A. Redshift Fluctuations 



In the average rest frame of the source and detector, let v M = t M = (1,0,0,0) and 
= (1,1,0,0) as previously. With this and the symmetry and cyclic properties of the 
Riemann tensor, Eq. (J73|) reduces to 

((AO 2 ) = J da J da 1 (R mm {x)R mm {x')) p . (76) 

We construct the thermal Riemann tensor two point function via the Matsubara method as 
an infinite image sum in imaginary time of the vacuum two point function. We proceed for 
the moment in the Feynman gauge, but since we are strictly dealing with functions of the 
Riemann tensor, the result will be gauge independent. From Eq. (|30|) and the symmetry 
properties of the metric tensor, the Riemann tensor vacuum two-point function can be 
written as 



(ifooi(z)-ftoioi(z')) 

= ^((2/?-oi,oi — — ^n,oo) ( x ) (2/ioi,oi — ^00,11 — ^-11,00) i x ')) 

= -[4Goioi,oioi — 2Goioo,oin — 2Goin,oioo — 2Goooi,noi + ^0000,1111 

+ GWl.llOO — 2^1101,0001 + ^iioo.OOll + Gim.oooo]- (77) 

Here, G a p^ v is the metric two point function. The first two indices on G Q ^ refer to point 
x while the second two refer to point x' and similarly for the derivative indices. In the 
Feynman gauge 

p%v)D (78) 

where D is the vacuum scalar two point function 

D = — - (79) 

4vr 2 (Ax) 2 - (At) 2 ' 1 J 

One finds that the only nonzero components of are G ioi, GW>0; GWi, Guoo, and Gnu 

and thus 

(-Roioi( a; )-Roioi(^ / )) = ~~^c~ ^D,txt'x' +D, XXX ' X / +D, xxt ' t / +D, ttx ' x ' +D, m >t> )■ (80) 

One can easily see from the function D that D )tx t'x> = D, xxt > t / = D, ttx ' x > = d^d^D and there- 
fore 

<itW(^oioi(*')> = \{dt - 2d 2 t d 2 x + &t)D. (81) 

It is now clear that the thermal Riemann tensor two point function may be constructed from 
the vacuum two point function by making the replacement D — > Dp whence 

(iW^oioi (*')>/? = \{St ~ 2d?d 2 x + dl)D p . (82) 
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As mentioned earlier, Dp is constructed via the Matsubara method. First make D periodic 
in imaginary time and then take an infinite image sum so that 



D 



oo 

—y 



1 



4vr 2 ^ (Ax ) 2 -(t-t' + in/3) 2 ' 

n=— oo 



(83) 



The prime on the summation indicates that we leave out the n = term, which corresponds 
to the vacuum term. As a result, Eq. (jH^j) gives a normal-ordered two-point function. We 
can set y — y' = z — z' = 0, and use null coordinates to write 



oo 

-y 



1 



(84) 



This, with the null coordinate version of Eq. (JS2J), gives 

((AO 2 ) =4 / da J dd d u d u ,d v d v ,D 



a- 



(85) 



The integral can be evaluated in the following way. First note that we are interested in the 
real part of Dp, which is even in n, so that we may make the replacement 



«J oo 



n=l 



and let 



where 



1 oo 



n=l 



G = Re 



1 



n 2 /3 2 - AvAu - in(3(Av + Au) 



In null coordinates, 



9 °° r^o 

n=l 1/0 



to 



du \ dv! 



o 



m+2s 



Performing the integrations on t> and v' yields 



dv 



u'+2s 



dv' d v d v 'd u d u 'G. 



(86) 



(87) 



(88) 



O ^ /-to /"to 

((A0 2 ) = -$^ / ciu / du'I(Au) 
k n=1 Jo Jo 

where we make use of the fact that d u d u /G = —d\G in writing 



I(Au) = (dlG) 



Au=Aw+2s 



+ (d 2 u G) 



Av=Au-2s 



2(9 2 G) 



Av=Au 



(89) 



(90) 



The function /(Au) is a function of Au only, for which 

/•to rto rUt i-o 

du du'I{Au)= d(Au)(t - Au)I(Au) + d{Au)(t + Au)I(Au). (91) 

JO Jo Jo J -t 
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This now leads to the expression 

((A0 2 ) = — YW d(Au)(t Q -Au)I(Au) + d(Au)(t + Au)I(Au). (92) 
^ n=1 ./o ./-to 

Equation (J9~2J) can be evaluated using a program such as Maple. The complete result is 
a rather lengthy expression and no insight is gained by writing it out. However, in the 
limit where s ^> t and s ^> (3, i.e., the observationally reasonable limit where the distance 
between source and detector is much larger than both the observation time and the thermal 
wavelength, the expression reduces to 

((A£) 2 ) ~ — 2 i— [STrVcschV) + ttV - 3], (93) 

with a = t //3. In the limit where the observation time is short compared to the correlation 
time P, i.e. t <C /3, we find 

U ^ ; 45/3 4 45/? 4 { ' 

In this case, the rms line width grows linearly with observation time, 



While on the other hand if t 3> /?, 



Here the rms line width approaches a constant, 



(AOrrns = ^ (97) 



B. Angular Blurring 

Turning our attention to Eq. (JZEJ, once again let = t M = (1, 0, 0, 0) and = (1, 1, 0, 0); 
additionally let = (0,0,1,0). With this and the symmetry and cyclic properties of the 
Riemann tensor, Eq. (|75|l reduces to 

((A0) 2 ) =jdaj da' ((i? 210 i(x) - R W ia(x)){R 2101 (x') - R 20 io(x')))p. (98) 

Proceeding in the Feynman Gauge, the calculations follow those for redshift fluctuations. It 
is straightforward to show that 

{{Ryxtx{x) ~ Rytxt{x)){Ryxtx{x') ~ Rytxt{x'))) 

= - A {d$ + 2dld x - 2d t dl - dl)D = -Ad u d 2 v d v ,D . (99) 
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In null coordinates, the equation for angular blurring becomes 



ptQ ptQ ru-\-Zs ru -\-Zs 

((AG) 2 ) = -4/ du du 1 dv dv'd u d 2 v d v ,D. (100) 

JO JO Ju Ju' 

This integral can be solved via the same method invoked in the previous subsection when 
dealing with redshift fluctuations. While it is not immediately obvious from this form of the 
integral, the limiting results are the same as for line broadening. Namely, in the case when 
s ^> t , p we again have 

<( A0 ) 2 > ~ 5-iW 3 * 2fl2csch2 M + ~ 3 ]' ( 101 ) 

c/7T (X ID 

In the limit where a — > 0, or to <C /3, 

(AeU.-^, (102) 



/3 2 V45 ' 

while on the other hand when a — > oo, or to ^ A 

Of 

(A0) rms ~ . (103) 

As emphasized earlier, these results are independent of gauge choice. As a check, these 
results were also calculated in the TT gauge using the thermal two point function and 
Hadamard function for the TT gauge derived in Ref. 1231 . 

We may compare Eq. ()103|) with the results of Ref. [l2|, where a heuristic result is found 
for the angular blurring of the image of a source caused by passive fluctuations rather than 
active fluctuations as done here. The source of fluctuation there is taken to be thermal 
fluctuations in the stress tensor of a scalar field. In the high temperature limit, the passive 
fluctuation result is 

(A6) pass ~ , (104) 
bp 2 

where b is a characteristic width of the bundle of rays. The effect is seen to increase with 
source-detector separation, though it should be mentioned that this result is valid for the 
regime s < f3 3 /£py- An important point is that the result for passive fluctuations has 
two powers of £p\, while those for the active fluctuations have only one. One may therefore 
conclude that the effect tends to be larger for active fluctuations than for passive fluctuations. 



V. SUMMARY AND DISCUSSION 

In this paper, two effects of spacetime geometry fluctuation are examined. Expressions 
are derived for fluctuations of redshift and angular position of a source as observed by 
a detector. The physical manifestation of these effects is found to be a broadening of 
observed spectral lines and angular blurring of the image of a source. The fluctuation 
of these observables depends on the Riemann tensor correlation function, which in turn 
characterizes fluctuations in the spacetime geometry. These effects should arise for both 
active and passive metric fluctuations. However, in the case of passive fluctuations, it may 



18 



be necessary to do some additional spacetime averaging, as was discussed in Ref. jl2|. The 
explicit examples discussed in the present paper concerned active fluctuations. 

The effects of active spacetime fluctuations are examined by considering a linearized 
model of quantum gravity with gravitons occupying a squeezed vacuum state. In the ab- 
sence of squeezing there is no effect, but the effects grow exponentially with the squeezing 
parameter and so in principle can be made quite large. The redshift and angular position 
fluctuations are related to the energy density of the wave, and finally an order of magnitude 
is estimated for a gravity wave with large energy density. The result of this estimation is 
quite small with questionable observability in the present day universe. The results, however, 
indicate this is an interesting model for investigating the properties of fluctuating spacetime 
geometries. 

Further, a thermal bath of gravitons is considered as the source of spacetime geometry 
fluctuation. An expression for spectral line broadening is provided for the limit that the 
source-detector separation is much larger than both the observation time and the thermal 
wavelength. For observation times short compared to the thermal wavelength, the rms 
spectral line width increases linearly with the observation time. On the other hand, for 
observation times long compared to the thermal wavelength, the rms line width approaches 
a constant. The results for angular blurring mirror the results for spectral line broadening 
in the same limit of large source-detector separation. 
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APPENDIX A: SQUEEZED STATES 

A squeezed vacuum state is the natural state for a quantum mechanically created particle 
occupying an in-vacuum state represented in an out-Fock space. Squeezed quantum states 
are generated via the unitary displacement and squeeze operators. Here we provide a brief 
summary of the relevant ideas and results for squeezed states following primarily the notation 
as found in Ref. 24]|; see also Refs. |25|, |2(| |23|. The displacement operator generates the 



set of coherent states and is defined by: 

D{a) = exp(aa t - a*a) = e aa \~ aa , (Al) 



1 2 



where a is an arbitrary complex number. We will here be primarily interested in generating 
the coherent state \a) by displacing the vacuum state |0): 

\a)=D(a)\0). (A2) 

The displacement operator transforms a and a' as: 

(a)aD(a) = a + a (A3a) 
D\a)a ] D(a) = a) + a* (A3b) 
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The squeeze operator is denned as 

5(C) = exp[icV-lc(a t ) 2 ], C = re", (A4) 

where ( is an arbitrary complex number. We use the convention in 24]| for the squeeze 
operator transformations of a and a): 

S' (C)aS(() = a cosh r — a)e ie smhr (A5a) 

S^(()a^S(() = a+coshr - ae^sinhr (A5b) 

The squeezed state \a, () is now obtained by squeezing the vacuum and then displacing it: 

\a,C) = D(a)S(C)\0). (A6) 
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